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INTRODUCTION 
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Methods  of  specifying  optical  extinction  coefficients,  k,  (visible 
through  infrared)  in  the  marine  boundary  layer  through  the  use  of  meteoro¬ 
logical  observables  are  needed  to  support  naval  electro-optical  systems.  The 
k-coeff icient  is  a  function  of  wavelength,  X,  and  atmospheric  water  content  in 
both  condensed  form  and  vapor.  The  concentration  of  condensation  in  the  form 
of  aerosols  is  usually  quantitatively  depicted  by  a  spectrum,  n(r),  where  r  is 
aerosol  radius.  Coefficient  k(X)  is  critically  dependent  on  the  distribution 
of  the  total  particle  number,  N,  within  n(r). 

Most  data  on  n(r)  in  the  marine  boundary  layer  have  been  obtained  near 
the  surface.  Published,  readily  usable  models  specifying  nfr)  (eg,  Wells  et 
al,  1977,  and  a  Navy  model  described  by  Noonkester,  1980)  have  been  developed 
primarily  from  che  surface  data.  A  recently  developed  model  specifying  n(r) 
and  k(A)  (Gathman,  1983)  has  been  accepted  as  a  replacement  of  an  earlier  Navy 
model  (described  hy  Noonkester,  1980)  for  naval  EG  applications.  This  Navy 
aerosol  model,  hereafter  called  "Navy  model"  or  "model,"  was  based  primarily 
on  near-surface  data,  is  readily  usable,  and  is  being  tested  for  various 
meteorological  conditions. 

Some  Nivy  aerosol  systems  will  depend  on  optical  paths  that  have  an 
appreciable  portion  above  the  surface;  and  to  be  acceptable,  a  model  must  be 
successful  for  the  above- sut  ..ace  needs.  Recent  acquisition  of  data  by  NOSC 
(Noonkester,  1982a,  1982b)  provides  an  opportunity  to  test  the  model  for 
above-surface  optical  paths.  These  data  are  detailed  observations  of  n(r)  at 
many  levels  in  maritime  stratus  cloud  layers.  Because  stratus  clouds  are 
common  over  large  regions  of  the  ocean,  capping  layers  containing  numerous 
aerosols,  they  pose  a  likely  limiting  environment  for  systems.  The  purpose  of 
this  document  is  to  compare  k(A)s  specified  by  the  model  with  k(A)s  calculated 
from  the  observed  nir)s  beneath  stratus  clouds.  A  comparison  of  n(r)s  speci¬ 
fied  by  the  model  with  observed  n(r)s  is  made  to  isolate  reasons  for  differ¬ 
ences  between  the  modeled  and  observed  k(A)s.  These  comparisons  indicate  that 
the  model  cannot  specify  nit)  end  k(A)  beneath  stratus  clouds.  Apparently  the 
model  must  be  modified  appreciably  before  it  can  be  applied  to  conditions 
beneath  stratus  clouds. 
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2.  NAVY  AEROSOL  MODEL 

The  Navy  model  {Gathman,  1983)  conjists  of  a  three-mode,  log-normal 
distribution  of  n(r)  used  in  conjunction  with  Mie-theory  calculations  to 
specify  k(i)  for  40  Is  in  the  range  0.20  um  <  X  <  40.0  pm.  The  real  and 
imaginary  parts  of  the  refractive  index  are  specified.  The  model  alio  speci¬ 
fies  the  water  vapor  absorption  for  the  40  Xs,  although  that  characteristic  is 
not  used  here. 

I.tt-  basic  contents  of  the  model  for  the  condensed  portion  of  the 
atmospheric  water  content  are  as  fellows: 


k(X)  =  10  ^  if  ^Q(X,  m)  n(r) 


r2  dr 


3 

Z  A, 

1  =  1  * 


.  2  /  r  \ 
~ln  Fr  ' 

l 


where  n(r)  is  in  cm  ^  pn  1  and  k(X)  is  in  km  1  .■  Also 


=  nondimensional  extinction  cross  section  for  spherical  particles 

=  refractive  index,  both  real  and  imaginary  parts 

2 

=  2000(AM)  [AM  =  air  mass  parameter,  in  values  of  1-10.  AM  =  1  for 

open  ocean;  AM  =  1 0  for  coastal  region.} 

=  the  maximum  of  0.5  or  5.866(w  -  2.2) 

-5 

=  the  maximum  of  1.4  x  10  or  0.01527(v'  -  2; 

=  24-hour  average  surface  wind  speed,  in  m/s 
=  current  surface  wind  speed,  in  m/s 


2  -  f 
6(1  -  f) 


(growth  factor) 


=  dry  particle  size  for  mode  i 
-2 

=  3.0  x  10  Mm 
-1 

=  2,4  x  10  urn 


=  2.0  Mm 

=  fractional  relative  humidity 

,  A^,  and  correspond  respectively  to  modes  1,  2,  and  3. 
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One  convenient  form  of  k(i,  t)  is  as  follows: 


k(X,  f) 


-3 

10  TT 

F  (  f  ) 


3  C. <  X,  f  ) 

101 


(3) 


Tabulations  of  Ch  are  available  for  the  40  Xs  at  fs  of  0.50,  0.85,  0.95,  and 
0.99.  Absorption  by  water  vapor  was  not  considered  here. 

When  w  and  w'  are  not  known,  the  model  provides  "default"  values  as 
fol lows : 


Tropical 

W 

=  w‘ 

=  4. 1  m/s 

Mid-latitude  summer 

w 

=  w’ 

=  4.1  m/s 

Mid-latitude  winter 

w 

=  w 1 

»  10.3  m/s 

Subarctic  summer 

w 

=  w 1 

=  6.7  m/s 

Subarctic  winter 

w 

=  w' 

=  12.4  m/s 

Although  the  Navy  model  was  designed  for  near-surface  conditions,  the 
model  was  applied  to  above-surface  stratus  layer  conditions  by  using  f  as  an 
independent  variable.  Because  f  errors  were  unknown  during  the  stratus  mea¬ 
surements,  the  f  measurements  were  not  used  to  establish  elevations  for  com¬ 
parison.  Although  an  adiabatic  vertical  profile  of  f  was  used  to  provide  the 
elevation  dependence  for  the  model,  this  assumption  is  not  critical  for  these 
comparisons . 


3.  STRATUS  LAYER  DATA 


a.  SENSORS 


Measurements  made  aboard  a  twin-engine  Piper  Navajo  flying  at  54  m/s 

included  elevation,  z,  temperature,  T,  dew  point,  T^,  and  n(r).  The  combined 

use  of  radar  and  pressure  altimeters  reduced  measurement  errors  in  z  to  less 

than  +3  m  up  to  700  m.  The  total  system  accuracy  of  the  T  and  T,  measurements 

d 

has  not  been  finalized.  A  PMS  ASSP-100  spectrometer  provided  n(r)  over  the 
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range  0.23  ym  £  r  £  14.7  pm,  and  a  PMS  OAP-200  spectrometer  provided  n(r)  over  fl 

the  range  14.2  pm  £  r  <  150  pm.  T  and  z  were  sampled  every  5  s  and  a  complete 
n(r)  was  obtained  every  8  s. 

b.  DATA  ACQUISITION 

Slow  descents  were  made  through  extensive  stratus  cloud  decks,  to  esti¬ 
mate  the  elevations  of  the  cloud  tops,  z^_  .  and  clc  .a  bases,  x^.  Horizontal 
runs  of  2  minutes  (6.44  km)  were  then  made  at  about  the  following  elevations: 

Zg  (near  surface) 


N 

04 

O 

c 

t  n 

c 

“c 

80 

m 

z 

c 

60 

m 

z  - 
c 

40 

m 

z 

c 

20 

m 

z 

c 

z  + 

c 

20 

m 

z  + 

40 

m 

(z  +  z  )/2  [midcloud] 
c  t 

z^  -  40  m 
t 

z 

t 

z  t  40  u 
t 

All  rms  were  made  both  into  and  with  the  estimated  surface  wind. 

The  maximum  and  minimum  standard  deviations  of  z  were  respectively  9.3 
and  1.6  m,  and  the  average  standard  deviation  was  4.7  m.  All  data  along  a 
horizontal  run  were  accepted  if  they  were  acquired  at  an  elevation  within  +7  m 
of  the  average. 
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SPECTRAL  DATA 


c. 


Measurements  of  the  number,  n,  of  aerosols  in  a  volume  of  1  cm  for  a 
bandwidth  of  1  pm  centered  at  a  radius  ware  made  every  8  s  during  each 
2-min  horizontal  run.  A  complete  spectrum,  n(r),  consisted  of  47  n(r^  )s.  A.n 
average  n(r^)  in  each  radius  band  was  obtained  by  averaging  all  n(r^)s  ob¬ 
served  during  each  2-min  run.  Each  8-s  spectrum  represented  429  m,  and  the 
2-min  average  n(r)  represented  6.44  km,  or  fifteen  b-s  spectra. 

Parameters  calculated  from  the  average  r.(r)  include  (1  )  the  total  number 
of  particles,  N,  (2)  the  mean  radius,  r,  (3)  the  liquid  water  coe.tent,  w* ,  and 
(4)  k ( A  >  for  Xs  of  0.53,  3.75,  and  10.59  pm.  The  following  equation  was  used 
to  calculate  w*: 


w« 


dr 


(4) 


where  D  is  the  density  of  pure  water,  in  g/’em  ,  k(X)  was  calculated  from 
equatiuu  (1).  The  aerosols  wore  assumed  to  be  spherical,  and  the  integration 
was  extended  over  the  range  0.23  pm  <_  r  _<  150  pm.  The  refractive  indices  used 
are  those  given  by  Selby  et  al  (1976).  The  model  used  another  set  of  refrac¬ 
tive  indices.  The  k(Xls  calculated  from  the  observed  n(r)s  will  be  called 
observed  k(X)s,  although  optical  extinctions  were  not  observed. 


Stratus  layer  aerosol  data  were  obtained  over  the  ocean  about  80  miles 
southwest  of  San  Diego  on  14,  28,  and  29  May  and  11,  13,  14,  17,  and  18  August 
1981.  Figure  1  gives  the  measurement  levels  and  the  cloud-top  and  cloud-base 
elevations  on  each  day  of  the  measurements.  The  cloud  base  v/3s  defined  to  be 
at  the  level  where  w*  =  0.02  g/m3,  as  discussed  in  section  5b. 


4.  AIR  MASS  FOR  STRATUS  DATA 


a.  N  AND  r 

Table  1  gives  N  and  r  at  the  surface  and  at  1 00  m  above  the  cloud  bases 
for  the  May  and  August  data.  N  was  greater  during  August  at  all  levels.  N 
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was  about  2.5  times  greater  during  August  at  1 00  m  above  cloud  base.  The  mean 
radius,  r(  was  larger  at  all  elevations  during  May.  Pruppacher  and  Klett 
(1978)  give  N  values  about  four  times  greater  in  continental  clouds  than  in 
marine  clouds  at  r  >  3  p  and  give  examples  showing  that  n(r)s  are  character¬ 
ised  by  large  rs  in  marine  air  masses.  Accordingly,  the  May  and  August  data 
could  be  considered  tc  represent  marine  and  continental  air  masses, 
respectively. 


Data 

Near  Surf 

aco 

1 00  m  above  Cloud  Base 

Ter ) od 

1  981 

N 

(cm  3  ) 

r 

(cm ) 

w* 

(g/in3) 

N 

(cm  3  ) 

r 

(cm) 

w* 

(g/m3 ) 

May: 

3  days 

62 

0.67 

-4 

3.6  x  10 

269 

4.0 

j 

0. 1  5 

August 

• 

:  5  days  j 

i 

84 

i 

0.  39 

9.4  x  10~5 

)  1 

665 

3.  1 

0.21 

Table  1.  Average  aerosol  spectrum  parameters  representing  a 
horizontal  distance  of  6.44  km. 


b.  MODE  IN  n ( r ) 

In  the  May  data,  distinct  modes  formed  in  n(r)  above  the  cloud  base.  The 
aerosol  radius  at  the  peak  of  the  mode  increased  linearly  with  elevation  from 
3  Mm  near  cloud  base  to  a  maximum  just  below  cloud  top.  Modes  did  not  form  in 
the  clouds  during  August,  apparently  because  a  much  greater  number  of  cloud 
condensation  nuclei  were  present  that  were  capable  of  competing  for  the  avail¬ 
able  moisture.  Nieburger  and  Chien  (I960),  Fitzgerald  (1974),  and  Lee  et  al 
(1980)  presented  cloud  models  that  produced  modes  in  n(r)  for  marine  air 
masses • 
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AIR  HASS  SOURCE 


C. 

The  measurement  region  was  near  118°W  longitude  and  32®N  latitude,  as 
shown  in  figure  2.  An  average  surface  pressure  map  was  constructed  for  the 
region  bounded  by  longitudes  from  110®  to  140®  west  and  latitudes  from  30®  to 
55®  north.  Pressures  were  taken  from  the  0400  PST  synoptic  maps  at  intersec¬ 
tions  of  5°  longitude  and  latitude  increments  within  this  region  for  the  days 
13,  14,  27,  28,  and  29  May,  to  represent  the  May  data,  and  10,  11,  12,  13,  14, 
16,  17,  and  18  August,  to  represent  the  August  data.  After  average  pressures 
were  computed  for  each  intersection,  maps  of  the  average  pressures  were  con¬ 
structed  for  the  May  and  August  days,  and  these  are  presented  in  figure  2. 

The  average  pressure  map  for  May  shows  an  offshore  subtropical  high  in  a 
southerly  position  and  a  low  adjacent  to  the  northwest  part  of  the  region. 
During  August  the  average  pressure  map  shows  a  subtropical  high  elongated 
north  to  south  in  the  western  portion  of  the  region,  producing  a  northerly 
flow  west  of  120°  longitude  at  all  latitudes. 

The  northerly  flow  at  all  latitudes  in  figure  2  for  tha  August  days 
provides  strong  evidence  that  a  continental  air  mass  would  be  present  during 
the  August  measurement  days.  Because  the  pressure  patterns  are  highly  similar 
for  the  May  and  August  days  along  the  coast  of  California,  a  continental  air 
mass  might  be  expected  at  the  measurement  site  during  May.  However,  the 
west-east  pressure  gradient  is  greater  during  May  west  of  San  Diego.  A 
stronger  cross-isobaric  flow  (more  westerly  wind)  is  expected  in  the  San  Diego 
region  when  the  west-east  pressure  gradient  increases.  The  stronger  west-east 
pressure  gradient  for  the  May  days  might  have  been  sufficient  to  maintain  a 
marine  air  mass  in  the  measurement  region  during  May. 

Taken  collectively,  the  above  data  give  strong  evidence  toward  the  con¬ 
clusion  that  the  aerosols  had  a  continental  source  during  the  August  days  and 
give  appreciable  evidence  toward  the  conclusion  that  the  aerosols  had  a  marine 
source  during  the  May  days.  Accordingly,  air  mass  parameters  (AM)  of  1  and  10 
will  be  assigned  to  the  May  and  August  days,  respectively.  Subsequent  dis¬ 
cussions  will  show  that  the  AM  value  assignments  were  not  critical  in  the 
comparisons . 
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Figure  2.  Representative  surface  pressure  patterns  for  measurement  days  in  May 
and  August  1981 .  The  measurement  aiea  location,  near  San  Diego,  is  shown.  The 
pressure  patterns  indicate  a  probable  continental  air  mass  prescnl  in  the  measure¬ 
ment  area  in  August,  a  probable  marine  ai:  mass  in  May. 


5.  PROFILE  OF  kU  )  BENEATH  STRATUS 

a.  RELATION  BETWEEN  k< X )  AND  w* 

When  k(X)  and  w*  are  determined  from  n(r)  by  means  of  equations  (1)  and 
(4)  respectively  arid  the  coefficients  a  and  b  in  the  expression 

k( X )  =  a(w*)b  <5: 

are  determined  by  statistical  regression  analysis,  the  coefficient  of  correla¬ 
tion,  p,  between  k(X)  and  a(w*)b  is  found  to  be  greater  than  0.95  for  a  great 
variety  of  spectral  shapes.  The  coefficients  a  and  b  vary  with  X  and  appear 
to  vary  with  air  ma3s  (eg,  Pinnick  et  al,  1979s  Hughes  and  Jensen,  1978; 
Noonkester,  1980). 

Comparisons  of  the  modeled  and  observed  k(X)s  were  made  for  Xs  of  0.53, 

3.75,  and  10.59  pm.  Table  2  gives  a,  b,  and  p  for  these  Xs  determined  from 

3 

n(r>  stratus  data  except  n(r)  where  N  was  less  than  iQ/cm  ,  with  the  May  and 
August  days  considered  separately.  The  large  ps  indicate  that  w*  can  almost 
exactly  specify  k(X)  in  equation  (5). 


Data 

Period 

X 

( pm ) 

a 

b 

Correlation 

Coefficient, 

P 

May:  3  days 

0.53 

1  94 

0.834 

0.99 

3.  75 

308 

0.950 

0.99 

10.59 

308 

1.16 

1  .00 

August:  5  days 

0.53 

231 

0.796 

0.99 

3.75 

351 

0.949 

0.98 

10.59 

282 

1.13 

1  .00 

Table  2.  Values  of  a  and  b  in  equation^ 5)  determined  by 
regression  analysis.  k(X)  has  units  km  for  these  constants 


b.  AVERAGE  PROFILE  OF  w*  AND  k(X)  BENEATH  STRATUS  BASE 

If  the  average  vertical  profile  of  w*  below  cloud  base  is  known  for  the 
May  and  August  data,  then  the  average  vertical  profile  of  k(X)  below  cloud 

base  can  be  obtained  from  equation  (5)  by  using  values  of  a  and  b  in  table  2. 

The  bases  of  the  stratus  clouds  could  not  be  uniquely  determined  visually 
because  the  decrease  in  horizontal  visibility  from  "good"  to  "poor"  was 

gradual  and  extended  over  depths  ranging  from  40  to  100  m  during  ascent  into 

the  clouds.  The  bases  cf  the  stratus  clouds  were  assumed  to  be  at  the  level 
where  w*  =  0.02  g/m3.  When  the  average  cross-sectional  aerosol  area  at  this 
level  and  Koschmieder ' s  equation  are  used,  the  level  chosen  for  cloud  base  has 
a  visibility  of  435  m.  According  to  the  international  visibility  code,  this 
visibility  corresponds  to  a  moderate  fog. 

3 

Figure  3  presents  the  average  w*  below  the  cloud  base  (w*  *  0.02  g/m  ) 
for  the  May  and  August  days.  Distances  relative  to  this  defined  cloud  base 
will  be  identified  by  z* . 


Figure  3.  Profile  of  the  average  liquid  water  content  for  3  stratus  days  in  May  198)  (marine  air  mass)  and 
for  5  stratus  days  in  August  1981  (continental  air  mass). 
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The  minimum  w*  has  been  reported  to  range  from  0.01  to  0.05  g/m  near  the 
base  of  stratiform  clouds  (Houghton,  1951)  and  in  fogs  (McCartney,  1976; 

Heaps,  1962).  Using  this  range  of  w*s,  the  cloud-base  heights  would  be  in  the 
range  -10m<z*<20min  figure  3.  The  increase  in  w*  above  z*  *  20  m  is 
linear,  commensurate  with  moist  adiabatic  cooling  and  cloud-top  entrainment. 
Apparently,  the  level  at  which  w*  =  0.02  g/m"^  is  near  the  saturation  level  (f 

3 

=  1.0).  The  level  where  w*  =  0.02  g/m  is  considered  to  be  within  +15  m  of 
the  true  saturation  level. 


In  figure  3,  the  larger  values  of  w*  below  cloud  base  for  May  would  be 
expected  because  r  is  greater  at  all  levels.  The  near-constant  value  of  w* 
between  z*  =  -20  m  and  z*  =  -40  m  for  May  was  caused  by  a  perturbation  in  the 
data  at  one  level  on  29  May  and  is  not  considered  to  be  representative;  the 
dashed  portion  is  a  more  likely  w*  profile  during  May  in  this  region.  The 
more  rapid  increase  of  w*  in  the  region  -40  m  <  z*  <  0  during  August  is  con¬ 
sidered  to  be  caused  by  the  large  uptake  of  water  vapor  by  a  greater  number  of 
small  aerosols  as  saturation  is  approached. 

In  equation  (5),  k(X,  z*)  can  be  determined  by  using  w*(z*)s  from  figure 
3.  Figure  4  gives  the  observed  k(z*)s  for  is  of  0.53,  3.75,  and  10.59  pm 
separately  for  the  May  and  August  days. 

6.  PROFILE  OF  RELATIVE  HUMIDITY 

The  model  is  effective  in  providing  k(X)  as  a  function  of  elevation,  z, 
if  the  fractional  relative  humidity,  f,  can  be  specified  as  a  function  of 
elevation.  Thus,  if  f  earn  be  specified  as  a  function  of  distance  below  z*, 
the  k(X)s  given  by  the  model  can  be  compared  with  the  k( X,  z*)s  given  in 
figure  4. 

Observations  and  theory  show  that  the  region  below  the  cloud  base  approx¬ 
imates  a  well-mixed  adiabatic  layer.  Because  errors  in  the  fs  measured  below 
the  stratus  clouds  along  the  runs  are  unknown  and  may  be  large  at  f  >  0.9,  the 
deviations  of  the  actual  f(z*)s  from  those  in  an  adiabatic  layer  cannot  be 
determined.  An  adiabatic  lapse  rate  for  f  will  be  assumed  to  be  present  below 
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figure  4.  Profiles  of  optical  extinction  coefficients  by  aerosols  for  observations  (solid)  and  for  the  Navy 
model  (dashed).  In  the  Navy  model,  profile  A  represents  open  ocean  (AM  =  1)  with  4.1  m/s  wind  speed, 
profile  0  represents  coastal  (AM  =  1 0)  with  4.1  m/s  wind  speed,  and  profile  C  represents  either  open 
ocean  (AM  -  1)  or  coastal  (AM  =  10)  with  10 3  tn/s  wind  speed. 
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cloud  base,  and  f  will  be  assumed  to  be  1.0  at  the  cloud  base.  The  z*s  for  fs 
of  0.99,  0.95,  and  0.85  in  an  adiabatic  layer,  corresponding  to  fs  for  the 
tabulated  C^s  in  equation  (3),  are  respectively  -12,  -65,  and  -215m.  Data 
shown  in  a  later  section  will  demonstrate  that  large  changes  in  the  profile  of 
f  would  not  permit  the  model  to  approximate  the  observed  k(A)s;  thus  an 
assumption  of  an  adiabatic  lapse  rate  for  f  is  reasonable. 

7.  SURFACE  WIND  SPEED 

The  24-hour  average  surface  wind  speed ,  w,  and  the  current  surface  wind 
speed,  w',  are  the  two  remaining  parameters  unspecified  in  the  model.  Obser¬ 
vations  by  the  author  during  the  low-level  horizontal  runs  indicated  that  w‘ 
ranged  from  about  5  to  12  knots  (2.5  to  6  m/s)  on  all  8  days.  The  synoptic 
patterns  near  the  observation  days  indicated  that  the  local  wind  speed  should 
be  representative  of  a  large  region;  thus  the  current  and  24-hour  average  wind 
speed  can  be  assumed  to  be  the  same.  Acceptance  of  the  default  wind  speeds  w 
and  w'  of  4.1  m/s  for  a  midlatitude  summer  appears  reasonable. 

8.  PROFILES  Or'  MODEL  k(A)s 

a.  GRAPHICAL  PRESENTATION 

Figure  4  presents  model -provided  k(A,  z*)s  appropriate  for  an  adiabatic 
lapse  rate  for  f  below  the  stratus  cloud  base.  The  model  k(A,  z* )  profiles 
are  presented  separately  for  Is  of  0.55,  3.75,  and  10.59  pm  and  for  the  May 
and  August  data.  For  comparison,  model  K(A,  z*)s  are  given  for  AMs  of  1  and 
10  on  all  figures  (although  AMs  of  1  and  10  are  considered  to  be  appropriate 
for  May  and  August,  respectively).  The  ks  calculated  from  the  stratus  data 
for  A  =  0.53  urn  should  be  comparable  to  ks  tabulated  for  the  model  at  A  =  0.55 
lira.  Hereinafter,  the  observed  k  for  A  =  0.53  pm  will  be  considered  applicable 
to  A  =  0.55  pm. 


. ;,w' 
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b.  VARIATION  OF  k(X)  WITH  AIR  MASS 

The  difference  between  the  model  ks  for  AMs  of  1  and  10  increases  with  an 
increase  of  f  and  with  a  decrease  in  X.  The  differences  are  appreciable  for  X 
=  0.55  pm.  The  differences  at  Xs  of  3.75  and  10.59  ym  may  be  insignificant 
for  many  purposes. 

c.  VERTICAL  GRADIENT  OF  MODELED  k(X) 

The  vertical  gradients  (Ak/Az*)  for  the  model  were  calculated  for  the 
change  in  k  between  fs  of  0.99  and  0.95  and  between  fs  of  0.95  and  0.85.  The 
vertical  gradients  were  greater  for  the  larger  fs  and  for  the  AM  of  10.  For 
either  AM,  the  gradients  decreased  in  the  following  order  of  Xs:  3.75,  0.55, 
10.59  ym.  The  variation  with  X  for  an  AM  of  1  was  minor.  The  greatest  in¬ 
crease  in  the  vertical  gradient  of  k  for  any  X  was  16;  it.  occurred  for  a  X  of 
0.55  ytn  when  AM  changed  from  1  to  10  at  the  large  fs. 

d.  RELATIVE  CONTRIBUTION  OF  EACH  MODE  TO  k(X) 


Comparison  of  the  relative  contribution  of  each  mode  of  the  aerosol  model 
to  the  extinction  coefficient  k(X)  permits  evaluation  of  the  significance  of 
the  air  mass  factor  through  the  use  of  A,) ,  the  average  24-hour  wind  speed 
through  the  use  of  A^ ,  and  the  current  wind  speed  through  the  use  of  . 

Table  3  gives  the  percent  contribution  of  each  mode  to  k(X)  for  Xs  of  0.55, 
3.75,  and  10.59  at  fs  of  0.99,  0.95,  0.85,  and  0.50  in  an  open  ocean  (AM  =  l) 
and  a  coastal  environment  (AM  ~  10)  when  the  surface  wind  speed  (current  and 
24-hour  average)  is  4.1  m/s. 

Data  in  table  3  indicate  the  following: 

(1)  For  AM  =  1 ,  mode  3  is  most  important  and  mode  1  may  be  ignored. 
For  AM  =  10,  mode  1  controls  k(X)  for  X  =  0.55  ym  and  all  three 
modes  must  be  used  for  Xs  of  3.75  and  10.59  ym. 
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A 

(  Pm) 

f 

Open 

: 

Ocean  (AM  -  1  ) 

Coastal  (AM  = 

10) 

Contribution  (%) 

Contribution  (%)* 

Mode 

1 

Mode 

2 

Mode 

3 

Mode 

1 

Mode 

2 

Mode 

3 

0.55 

0.99 

1  3 

36 

51 

(94) 

(3) 

(4) 

0.95 

9 

39 

52 

(91  ) 

(4) 

(5) 

0.85 

6 

41 

53 

(86) 

(6) 

(8) 

0.50 

4 

43 

53 

(81  ) 

(9) 

(10) 

3.75 

0.99 

1 

43 

56 

(47) 

(23) 

(30) 

0.95 

0 

35 

65 

21  i 

28 

51 

0.85 

0 

27 

73 

ii  ! 

24 

65 

0.50 

0 

20 

80 

6 

19 

75 

10.59 

0.99 

0 

1  7 

83 

(31  ) 

(11  ) 

(58) 

0.95 

0 

.0 

90 

21 

8 

71 

0.85 

0 

8 

92 

16 

6 

78 

0.50 

0 

5 

95 

9  1 

5 

86 

•Parentheses  are  used  to  indicate  where  mode  1  contributes  more  than  25%  to 
total  k. 
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Table  3,  Percent  contribution  of  each  mode  to  k{\)  for  the 
Navy  aerosol  model  in  open-ocean  and  coastal  environments 
for  a  wind  speed  of  4.1  ra/s  (aerosols  only). 

(2)  The  change  in  the  contribution  of  the  modes  as  AM,  f,  and  X 
change  is  less  distinct  for  mode  2.  For  either  AM,  mode  2 
increases  in  importance  as  f  decreases  at  A  =  0.55  urn  and 
decreases  in  importance  as  f  decreases  for  Xs  =  3.75  and  10.59 
pm.  Mode  2  contributes  about  the  same  for  any  AM  when  X  =  10.59 
pm. 

(3)  For  either  AM,  mode  3  increases  in  importance  as  f  decreases  and 
as  X  increases.  Mode  3  is  more  important  than  mode  2  for  any  f. 

In  general,  these  features  indicate  that  the  importance  of  the  smallest  aero¬ 
sols  (mode  1 )  increases  as  f  increases  and  as  X  decreases,  and  that  the  impor¬ 
tance  of  large  aerosols  (mode  3)  increases  as  f  decreases  and  as  X  increases. 
These  general  relations  are  expected.  The  relative  intermediate  importance  of 
mode  2  for  an  AM  of  10  may  not  have  been  expected. 
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Because  mode  1  Includes  many  aerosols  in  the  radius  region  below  the 
minimum  observed  radius  (r  •  0,23  ym),  comparison  of  the  observed  and  modeled 
k(X)  is  not  appropriate  when  the  relative  contribution  of  mode  1  to  k(X)  is 
appreciable.  An  appreciable  contribution  for  mode  1  is  arbitrarily  considered 
to  occur  when  the  relative  contribution  of  mode  1  to  k  is  25  percent  or  more 
of  the  total  k.  All  affected  modes  in  this  category  are  enclosed  by  paren¬ 
theses  in  table  3.  All  fs  and  Xs  in  table  3  are  acceptable  for  comparison  for 
an  open-ocean  environment.  In  a  coastal  environment,  only  Xs  of  3.75  and 
10.59  ym  at  fs  of  0.95,  0.85,  and  0.50  are  acceptable  for  comparison.  In 
figure  4,  the  data  points  of  the  modeled  k(X)s  in  which  mode  1  contributes  25 
percent  or  more  to  the  total  value  are  enclosed  by  squares  and  will  not  be 
considered  in  the  comparisons. 

9.  CRITERIA  FOk  ACCEPTANCE  OF  MODEL 

a.  "ACCEPTABLE"  TEST 

A  rigid  test  of  the  model  would  require  accurate  observations  of  n(r) 
along  with  accurate  measurements  of  f,  w,  and  w'  and  an  appropriate  method  of 
obtaining  AM.  Reliable  differences  between  the  observed  and  the  modeled  n(r)s 
or  k(X)s  could  then  be  determined.  The  significance  of  the  difference  between 
the  modeled  and  observed  k(X)s  could  be  estimated  relative  to  an  optical 
system  application  if  criteria  for  successful  performance  of  the  system  are 
quantified.  Application  of  this  procedure  is  unlikely  since  accuracies  of  the 
observed  n(r)s  are  generally  unknown,  a  method  of  obtaining  AM  is  not  defined, 
and  system  performance  limitations  by  errors  in  k(X)  are  not  usually 
available. 

b.  LIKELf  ERROR  IN  OBSERVED  k(X)s 

Simultaneous  measurements  by  colocated  PMS  aerosol  measuring  devices  at 
the  NW  tip  of  San  Nicolas  Island  during  a  5-day  period  iri  May  1979  when  atmo¬ 
spheric  conditions  varied  considerably  showed  that  the  k(X)s  calculated  from 
the  observed  n(r)s  differed  by  a  factor  of  2  to  3.  k(X)s  determined  by  a 
nearby  nephelometer  were  within  a  factor  of  2  to  3  of  the  k(X)s  calculated 


from  the  observed  n(r)s  (Jensen,  1980).  Thus  the  k(.\)s  in  figure  4,  obtained 
from  the  average  ws,  should  be  within  a  factor  or  2  to  3  of  values  expected 
from  measurements  by  similar  devices  and  could  be  within  a  factor  of  2  to  3  of 
the  real  values  in  a  representative  stratus-cloud  layer.  Accordingly,  the 
modeled  and  observed  k(X)s  can  be  considered  to  have  appreciable  differences 
when  they  differ  by  a  factor  of  more  than  3. 

c.  VERTICAL  PROFILES  OF  OBSERVED  k(\)s 

No  evidence  has  been  found  to  indicate  that  any  systematic  errors  in  the 
measurements  of  n(r)  by  the  same  PMS  spectrometer  vary  within  a  period  of 
several  weeks.  Changes  in  the  observed  n(r)  from  the  data  period  in  May  to 
August  can  be  considered  real.  The  vertical  profiles  of  the  observed  k(X)s  in 
figure  4  are  considered  to  approximate  the  true  profile  closely,  although  the 
magnitude  of  the  k(X)s  may  be  systematically  too  large  or  too  small  by  a 
factor  of  2  to  3. 

d.  MODEL  ACCEPTANCE  TEST  USING  PROFILES  OF  k(X)s 

The  true  profile  of  f  may  not  approximate  an  adiabatic  layer  (for  which  f 
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-  1  at  the  level  where  w*  =  0.02  g/m  ),  as  assumed  in  the  construction  of  the 
profile  of  the  modeled  k(x)s  in  figure  4.  The  modeled  k(X)  profiles  would  be 
acceptable  if  reasonable  changes  in  the  assumed  vertical  profile  of  f  would 
cause  the  modeled  k(X)s  to  be  within  a  factor  of  3  of  the  observed  k(x)  pro¬ 
files  at  all  levels  and  \s  and  for  reasonable  values  of  w,  w' ,  and  AM.  A 
vertical  gradient  oE  f  greater  than  adiabatic  wculd  not  be  reasonable,  whereas 
a  profile  of  f  less  than  adiabatic  might  be  acceptable  when  f  =  1  is  within 
+  1 5  m  of  z*  =■  0.  The  vertical  adiabatic  gradient  of  f  is  about  +1.4%/20  m  for 
the  temperatures  on  the  days  of  measurement. 

10.  COMPARISON  OF  MODELED  AND  OBSERVED  k(X)  PROFILES 

Appendix  A  is  a  detailed  comparison  stating  that  the  model  k(X)s  cannot 
simultaneously  approximate  (within  a  factor  of  3)  the  observed  k(x)s  for  all 
Xs,  both  months,  and  all  elevations  —  even  with  unreasonable  changes  m  the 


profile  Of  relative  humidity  and  surface  wind  speeds  (w  and  w' ) .  This  conclu¬ 
sion  applies  as  well  to  the  surface,  where  the  model  would  be  expected  to  be 
reliable. 

11.  COMPARISON  OF  MODELED  AND  OBSERVED  n(r) 

Comparison  of  the  modeled  and  observed  n(r)  for  selected  values  of  AM  and 
f  should  reveal  some  reasons  for  the  large  differences  in  the  modeled  and 
observed  k(x)s.  The  default  wind  speed  of  4.1  m/s  given  by  the  model  for  a 
midlatitude  summer  appears  reasonable  for  the  synoptic  patterns  prevalent 
during  May  and  August  1981.  Reasons  given  in  section  4  support  selection  of 
the  model  "open  ocean"  (AM  =  1 )  to  represent  the  May  data  and  "coastal"  (AM  - 
10)  to  represent  the  August  data. 

Figures  5  and  6  present  the  modeled  and  observed  n(r)  for  fs  of  0.99, 
0.95,  and  0.85.  Figure  5  is  for  the  May  data  and  figure  6  is  for  the  August 
data.  The  left  part  of  figures  5  and  6  provides  the  contribution  of  each  mode 

of  the  model  to  n(r)  at  an  f  of  0.99  for  the  May  and  August  data, 

respectively. 

Compared  to  the  observed  n(r),  the  modeled  n(r)  provides  considerably 
fewer  (as  few  as  1/10  the  total)  aerosols  roughly  in  the  region  0.4  pm  <  r  < 

10  pm  for  May  and  fewer  aerosols  (as  few  as  1/5  the  total)  roughly  in  the 

region  1  pm  <  r  <  10  pm  for  August.  The  model  overestimates  n(r)  for  r  >  20 

pm  at  f  =  0,99  for  both  months.  The  large  positive  departure  of  the  observed 
n(r)  above  the  model  for  May  at  f  «=  0.85  (surface)  in  the  region  r  >  40  urn  is 
appreciable . 

The  most  significant  and  consistent  differences  between  the  modeled  and 
the  observed  n(r)  are  roughly  in  the  radius  region  centered  between  1  and  3  pm 
for  May  and  between  2  and  4  pm  for  August.  The  observed  n(r)e  are  similar  for 
May  and  August  when  f  «=  0.95  in  the  region  r  >  3  pm.  n(r)  is  about  40 
cm  "*pm  1  at  r  “  3  pm  for  both  months.  Given  AM  “  10,  w  =  10.3  m/s  (20  knots), 
w'  —  15.5  m/s  (30  knots),  and  f  =  0.99,  the  model  yields  n(r)  (fer  r  =  3  pm)  = 
cm  ^pm  Thus,  the  model  appears  incapable  of  producing  the  observed  n(r)  in 


contribution  of  each  mode  of  the  mode!  to  nfr)  at  an  f  of  0.99  Is  shown  at 


the  region  near  r  =  3  ym.  For  the  same  values  of  AM,  w,  w' ,  and  f,  the  mod¬ 
eled  n(r)  approximates  the  observed  n(r)  at  r  =  7  ym  and  greatly  overestimates 
the  observed  n(r)  at  r  =  30  Pm.  Apparently  an  increase  in  A.]  and  A^  is  re¬ 
quired  to  bring  the  modeled  and  observed  values  to  about  the  same  value  in 
these  radius  bands. 

The  observed  n(r)  had  a  mode  in  n(r)  at  small  rs  for  large  fs;  whereas 
the  model  does  not,  unless  A1  can  be  assumed  to  vanish.  This  n(r)  modal 
characteristic  has  been  observed  by  others  (JG  Hudson,  University  of  Nevada, 
and  EE  Hindman,  Colorado  State  University)  and  can  be  assumed  to  be  real. 
Although  the  observations  cannot  provide  information  on  n(r)  for  r  <  0.23  pm, 
the  mode  in  n(r)  at  small  r  and  large  f  must  be  considered  as  an  error  source 
in  the  model  if  applied  for  above-surface  marine-stratus  conditions. 


The  modeled  and  observed  k(X)s  were  almost  identical  for  both  months  at  X 
=  0.55  ym  for  AM  =  1 0  and  a  wind  speed  of  4.1  m/s,  although  the  im  ’el  produces 
many  more  aerosols  at  small  rs.  Apparently,  in  the  determination  of  k  lor  X  = 
0.53  ym,  the  larger  observed  n(r)  in  the  region  0.4  yn  <_  r  <_  10  ym  compensated 
for  the  absence  of  small  aerosols  given  by  the  model  at  much  smaller  r. 

12.  SUMMARY  AND  CONCLUSIONS 

An  analytical  model  recently  formulated  by  the  Navy  (Gathman,  1983) 
specifies  the  extinction  of  optical  radiation  by  water  vapor  and  aerosols 
(equations  1,  2,  and  3),  for  surface  conditions.  The  controlling  parameters 
are  a  24-hour  average  surface  wind  speed,  w,  a  current  surface  wind  speed,  w' , 
an  air  mass  factor,  AM,  and  the  relative  humidity,  f.  The  model  was  developed 
from  somewhat  limited  surface  data  and  is  being  refined  as  new  data  became 
available. 

A  recently  acquired  unique  set  of  aerosol  data  beneath  mari?ie  stratus 
clouds  (Noonkester,  1982a,  1982b)  was  used  to  determine  the  capability  of  the 
model  to  specify  the  optical  extinction  coefficients,  k,  by  aerosols  beneath 
marine  stratus  (water  vapor  not  included).  The  applicable  elevations  of  the 
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model  were  established  by  assuming  various  vertical  profiles  of  f,  initially 
assumed  to  be  adiabatic.  Surface  wind  speeds  w  and  w'  given  by  the  model 
(default  value)  were  accepted  as  representing  the  stratus  days.  The  stratus 
data  were  divided  into  two  distinct  air  masses  —  marine  (AM  =  1)  and  conti¬ 
nental  (AM  =  10)  —  on  the  basis  of  the  average  pressure  patterns  representing 
the  days  of  measurements  and  aerosol  spectra  characteristics. 

•The  model  ks  were  compared  with  the  ks  calculated  from  the  observed 
spectra  for  elevations  where  f  is  0.99,  0.95,  and  0.85;  for  w  =  w'  =4.1  m/s; 
and  for  Xs  of  0.55,  3.75,  and  10.59  vim.  Differences  between  the  observed  and 
model  profiles  of  k  were  generally  excessive  for  a  near-adiabatic  profile  of  f 
and  the  above  values  of  w,  w1 ,  and  AM.  Differences  remained  excessive  where 
unlikely  profiles  of  f ,  unlikely  values  of  w  and  w' ,  or  intermediate  values  of 
AM  were  used.  No  combination  of  unlikely  profiles  of  f,  values  of  w  and  w' , 
or  values  of  AM  could  be  found  that  reduced  the  differences  between  the  ob¬ 
served  and  model  profiles  of  k  to  acceptable  magnitudes  (less  than  a  factor  of 
three)  for  most  of  the  profiles  of  k.  On  the  basis  of  these  comparisons,  it 
is  concluded  that  without  modification,  the  model  is  incapable  of  reproducing 
the  observed  ks  produced  by  aerosols  beneath  marine  stratus  clouds. 

Comparison  of  model  and  observed  aerosol  spectra,  n(r),  (r  is  radius) 
revealed  appreciable  differences  in  many  regions  of  r,  assuming  an  adiabatic 
profile  of  f,  w  =  w'  =4.1  m/s,  and  AMs  of  1  and  10.  Reasonable  changes  in 
the  profile  of  f,  in  w  or  w' ,  and  in  AM  appeared  incapable  of  significantly 
reducing  the  differences  in  all  regions  of  r.  The  largest  difference  was  in 
the  region  I  pm  <  r  <  10  pm,  where  the  observed  n(r)  had  many  more  aerosols  at 
all  fs,  particularly  for  an  AM  of  1  (marine  air  mass).  Without  modifications, 
the  model  appears  incapable  of  reproducing  the  observed  aerosol  spectra. 

Similar  comparisons  should  be  made  by  analyzing  stratus  cloud  data  having 
reliable  measurements  of  f,  w,  and  w' .  Methods  are  needed  for  establishing 
values  of  AM  easily. 
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APPENDIX  A:  COMPARISON  OF  MODELED  AND  OBSERVED  k(X)  PROFILES 


DIFFERENCES  WITHOUT  ADJUSTMENTS  OF  f (z*)  OR  CHANGES  IN  w  OR  w' 

Generally,  the  modeled  k(A)s  differ  significantly  from  the  observed.  The 
modeled  k(A)s  for  an  AM  of  1  differ  more  than  those  for  an  AM  of  10,  and  the 
difference  increases  with  an  increase  in  elevation  (greater  fs).  Nevertn. 
less,  k  (A  =  10.59  pm)  at  f  =  0.85  is  within  a  factor  of  2  from  the  observed, 
for  both  AMs  and  both  months . 

k (X)  AT  THE  SURFACE 

The  average  observed  f  at  the  surface  was  0.83  (at  31  m)  for  May  and  0.81 
(at  41  m)  for  August.  These  fs  are  about  0.09  larger  than  expected  in  an 
adiabatic  layer,  assuming  f  =  1  at  the  cloud  base.  The  modeled  k(A)s  for  an  f 
of  C.85  should  represent  the  observed  k(A)s  at  the  surface,  since  the  modeled 
k(A)s  are  somewhat  insensitive  to  f  for  fs  near  0.80.  Accordingly,  the  ob¬ 
served  and  modeled  k(A)  (f  =  0.85)  at  A  =  3.75  pm  for  August  and  at  A  =  10.59 

pm  for  May  were  essentially  equal.  All  other  modeled  k(A)s  at  f  =  0.85  dif¬ 

fered  from  the  observed  at  the  surface  by  factors  greater  than  3.  Thus,  the 
model  does  not  appear  capable  of  specifying  the  k(A)s  at  the  surface;  adjust¬ 
ments  of  w  and  w'  cannot  cause  the  modeled  k(X)s  to  approach  the  observed 

within  a  factor  of  3  for  all  as  and  both  months  simultaneously. 

The  model  specifies  larger  k(A)s  at  the  surface  during  August,  whereas 
the  observed  profile  indicates  the  opposite  relation.  As  shown  below,  no 
reasonable  change  in  the  model  parameters  will  produce  the  observed  profile. 

COMPARISON  WITH  CHANGES  IN  f(z) 

No  change  in  f(z)  could  make  the  modeled  k  approach  the  observed  either 
at  a  A  of  0.55  (considering  only  AM  =  1 )  for  both  months  or  at  a  A  of  3.75  pm 
for  May.  If  the  elevation  at  which  f  -  1  were  decreased  by  50  to  80  m,  the 
modeled  profile  of  k(A)  would  approximate  the  observed  for  A  =  3.75  pm  during 
August  and  for  A  =  10.59  pm  during  May.  Such  a  change  in  the  elevation  of 
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saturation  is  not  reasonable.  No  reasonable  change  in  f(z)  would  cause  the 
modeled  k(X)  to  approximate  the  observed  at  all  levels,  for  all  Xs,  and  for 
both  months  when  the  wind  speed  is  4.1  m/s,  simultaneously. 

COMPARISON  WITH  LARGER  w  AND  w' 

If  both  w  and  w'  are  10.3  m/s  (20  knots)  and  AM  =  1,  curve  C  for  X  = 

0.55  ym  (in  main  text  figure  4)  is  generated.  This  profile  of  k(A)  differs 
from  the  observed  by  more  than  a  factor  of  3  at  fs  of  0-99  and  0.95  and  by  less 
than  a  factor  of  3  at  f  of  0.85.  For  Xs  of  3.75  and  10.59  Vim,  curve  C  is  for  w 
and  w*  values  of  10.3  m/s  and  an  AM  of  10.  These  curves  approximate  the 
observed  at  all  fs  for  X  =  3i75  ym  in  August,  approximate  the  observed  for  X  = 
3.75  Jim  at  f  =  0.85  for  May,  approximate  the  observed  for  X  =  10.59  urn  at  f  * 
0.95  for  May,  and  approximate  the  observed  for  X  =  10.59  jim  at  f  =  0.99  for 
August.  The  difference  is  excessive  for  fs  <0.88  at  a  X  of  10.59  vm  for  both 
months.  Thus,  a  change  of  wind  speed  will  not  cause  the  modeled  k(X)  to 
approximate  the  observed  at  all  levels,  for  all  Xs,  and  for  both  months, 
simultaneously,  when  the  humidity  profile  is  adiabatic  below  the  cloud  base. 

COMPARISONS  WITH  ADJUSTMENTS  IN  f(z*),  w,  AND  w' 

If  the  elevation  at  which  f  =  1  were  decreased  by  50  to  70  m,  curve  C  in 
figure  4  would  approximate  the  observed  k(X)  at  X  ■=  0.55  jim  for  both  months 
and  at  X  =  3.75  jim  for  May.  Such  a  change  for  the  remaining  three  X-month 
combinations  would  increase  to  excessive  magnitudes  the  differences  between 
the  modeled  and  observed  profiles.  No  combination  of  simultaneous  adjustments 
in  f(z*),  w,  and  w'  will  cause  the  modeled  k(X)s  to  approximate  the  observed 
k(X)s  at  all  Xs  and  for  both  months. 
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